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We propose a new method to model cluster scaling relations in modified gravity. Using a suite of non- 
radiative hydrodynamical simulations, we show that the scaling relations of cumulative gas quantities, such as 
the Sunyaev Zel’dovich effect (Compton-y parameter) and the x-ray Compton-Y parameter, can be accurately 
predicted using the known results in the ACDM model with a precision of ~ 3%. This method provides a 
reliable way to analyze the gas physics in modified gravity using the less-demanding and much more efficient 
pure cold dark matter simulations. Our results therefore have important theoretical and practical implications in 
constraining gravity using cluster surveys. 


I. INTRODUCTION 

Clusters of galaxies are the largest gravitationally-bound 
objects in the Universe, and their formation and evolution are 
strongly affected by gravity. Therefore, their abundance and 
distribution are sensitive to the nature of gravity, making them 
an ideal cosmological probe to test and constrain gravity theo¬ 
ries m. With them being observed in various ways, including 
x-ray, optical richness, the Sunyaev Zel’dovich (SZ) effect Q 
and weak lensing, we are entering a new era of precision clus¬ 
ter cosmology. The upcoming eROSITA project Q, for ex¬ 
ample, will probe up to ~ 10® galaxy clusters on the whole 
sky out to redshift z > 1, and there will be follow-ups and 
synergies with other surveys, such as Euclid |l4), which will 
further improve the determination of cluster properties. 

However, the robustness of the cosmological constraints 
from cluster surveys depends strongly on the accuracy of the 
measurement of cluster mass. The latter is highly non-trivial 
and, indeed, is one of the main challenges in cluster cosmol¬ 
ogy. Although the mass of a cluster can be directly derived 
from the gas density and temperature profiles in x-ray sur¬ 
veys, this method requires high-quality spectra and therefore 
long exposure time, which is expensive for distant clusters 
(z > 0.5) and also for a large number of clusters. It also 
assumes hydrostatic equilibrium in clusters, which induces a 
systematic error of ~ 10% in the cluster mass estimation Q. 

In practice, this difficulty can be overcome by using scaling 
relations which relate the cluster mass to mass proxies that are 
combinations of observables which are easy to measure and 
have small intrinsic scatters. These scaling relations can be 
calibrated either by observations or by hydrodynamical simu¬ 
lations 0121 . In observations, this can be done by using either 
complementary observables, e.g., weak lensing, or a subset 
of observational data with better mass determination. In the 
ACDM model, with ever-improving resolution and modelling 
of baryonic physics, simulation calibration of the scaling re¬ 
lations are becoming more accurate and found to have a good 
match with observations 0. 

However, using cluster observations to test gravity in a fully 
self-consistent way is nontrivial. People usually compare the 
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observational data with theoretical predictions for certain rep¬ 
resentative modified gravity models 0. One important point, 
which is sometimes missed, is that the calibration of the ob¬ 
servational data itself in such models are generally much more 
nontrivial. The intrinsic non-linearity in these models not 
only changes the various scaling relations from the known 
results in ACDM, but also their intrinsic scatters and cor¬ 
relations with other cluster properties. For example, in the 
widely-studied /(i?) gravity model ||9l, the lensing and hy¬ 
drostatic masses of a cluster can be different, and the differ¬ 
ence depends on the screening effect which further depends 
on many factors such as the cluster’s redshift, mass, and envi- 
ronment(see Ref. Qoi for detailed discussion). This makes it 
difficult, both in observations and simulations, to calibrate the 
scaling relations that involve hydrostatic observables, such as 
the x-ray luminosity and temperature. 


To tackle this important but so far not well-understood is¬ 
sue, in this work, we propose a new method to model cluster 
scaling relations in modified gravity. We show that the ef¬ 
fect of modified gravity can be modelled as a rescaling of the 
cluster gas mass fraction. This rescaling can be quantified 
directly using less demanding pure dark matter simulations, 
and then used to accurately predict cluster scaling relations in 
modified gravity, based on the known results in ACDM. We 
demonstrate the accuracy of this method by comparing the 
predictions from our rescaling method for three mass proxies 
- the x-ray luminosity Lx, the SZ Compton-y parameter Fsz 
and the x-ray Compton y parameter lx im -with those from 
non-radiative hydrodynamical simulations. We find that the 
method works a ~ 3% accuracy for the latter two proxies. 


This paper is organised as follows: In Sec.[^ we introduce 


the hydrodynamical simulations used in this work. In Sec. HI 


we introduce the halo catalogs constructed in this work. In 
Sec. |IV[ we describe the features of profiles of the gaseous 
halos in our simulations. In Sec.jV] we compare the gas frac¬ 
tions in ACDM and f{R) gravity. In Sec. VI we show that 
the effect of modified gravity can be modelled as a rescaling 
of the cluster gas mass fraction. In Sec. |VII| we summarise 
and conclude this work. 
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II. HYDRODYNAMICAL SIMULATIONS 

We ran a suite of hydrodynamical simulations using ECOS- 
MOG 02 , a modified gravity simulation code based on the 
publicly available RAMSES ifTSll code. We simulated an f{R) 
model which exactly reproduces the ACDM background ex¬ 
pansion history ca, using a box of size Lbox = lOO/i ^Mpc 
that contains N = 256^ cold dark matter particles. The cos¬ 
mological parameters are the same as the Planck ifTSi best-fit 
ACDM model (0° = 0.045,0° = 0.271,0° = 0.684,/i = 
0.671, ris = 0.962, and erg = 0.834). Initial conditions were 
generated using the MPGRAFIC package lfT6l at z = 49. 

Due to the expensive cost of hydrodynamic simulations, we 
focus only on the f{R) model with fjio = —10“® (where 
fat) is the present value of d//di?). We run three realizations 
in total for the /(i?) model, and for each f{R) simulation, 
we run two ACDM simulations as control, one with exactly 
the same initial conditions as the corresponding f{R) run and 
the other with a rescaled baryon fraction, |n°, while keeping 
= r 2 ° -I- n° unchanged. 

We assume that all the baryon component in our sim¬ 
ulation is non-radiative ideal gas which obeys Pgas = 
^BTgasPgas/(M''np), where Pgas is the thermal pressure of 
gas, Pgas = pm^ngas is the density, ks the Boltzmann con¬ 
stant, fi the mean molecular weight, rigas the number density 
of gas particles and nip the proton mass. We also assume 
that the gas is fully ionised with compositions of electrons, 
hydrogen and helium. We take the primordial mass fraction 
of hydrogen as nH/(nH + 4nHe) ~ 0.75, so that the mean 
molecular weight is p, « 0.59 . Finally, for a monatomic gas 
the adiabatic index is 7 = |. 


III. GAS TRACER PARTICLES AND EFFECTIVE HALO 
CATALOGS 

RAMSES is a mesh-based code which does not define gas 
particles by default. In order to trace the gas density field on 
the simulation grids, we sample it with tracer particles fol¬ 
lowing a Poisson process. The tracer particles have a uniform 
mass and their local mean is set to be proportional to the lo¬ 
cal gas density on the grid. By comparing the density power 
spectra of the tracer particles and of the density field on the 
grids, we find that with a large enough number of particles 
(e.g. A"gas = 800° as shown in Fig. [U, the tracer particles 
can accurately represent the original gas density field on the 
simulation grids. 

Our halo catalogs are constructed using a modified version 
of the publicly available Amiga Halo Finder (Ahe) iflTl . We 
call the halo catalog produced using the true density field, the 
standard catalog. On the other hand, as discussed in ifTSl , the 
modified Poisson equation in f{R) gravity can be cast into the 
standard form 

= 47rGa2jpeff , (1) 

by defining the effective energy density Jpeff which incorpo¬ 
rates all the modified gravity effect, with G being Newton’s 
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FIG. 1. A comparison of the power spectra from the gas density 
field on AMR grids (black solid line) and the gas tracer particles 
with two different numbers of particles (stars). With a large enough 
number of particles (e.g. Wgas = 800° ), the tracer particles and the 
original gas density field have almost the same power spectra, which 
indicates that the tracer particles can accurately represent the original 
gas density field on the AMR grids. 


constant. The halo catalog constructed using JpeS is referred 
to as effective catalog (see for technical details). 


IV. PROFILES 

We use halos in the mass range of 10^^Mq/Ii < M < 
lO°° ®M 0 //i, in which most halos are unscreened for /j^o = 
— 10“°. This range is further divided into two different bins: 
10°°Mo//i < M < 10°°'%/o//r and < M < 

IQts.sM q/H. We consider profiles only within lOkpc/ft, < 
r < Rvir, where the virial radius i?vir is defined by an over¬ 
density of 300 with respect to the critical density. 

The gas temperature profiles are shown in the first row of 
Fig. i We can see that the effective halos in f{R) gravity 
and the ACDM halos are very similar to each other. This is as 
expected since gas temperature is closely related to the gravi¬ 
tational potential of halos. The physical reason behind this is 
that, during gas accretion, the shocked gas converts the energy 
that it gains from gravitational infall into the thermal energy 
of itself. Since the bulk motion of gas accounts for only a 
small fraction of the total gained energy, the gravitational po¬ 
tential energy is mainly converted into the its thermal energy, 
and so the temperature of the gas is mainly determined by the 
gravitational potential. 

Compared with the temperature profile, the density profile 
is more complicated. Unlike cold dark matter, gas has a core 
at the halo center and its density profile is flat at r < i?vir/20. 
The specific entropy of gas does not vary significantly in this 
core region, as shown in the third row of Fig. and the 
shocked gas can be considered as adiabatic 1^ . From sim- 
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FIG. 2. Profiles for different gas quantities from four different halo 
catalogs: ACDM (black lines), ACDM with rescaled baryon fraction 
I fit. (green lines), standard halo catalog in f{R) gravity (blue lines), 
and effective halo catalog in f{R) gravity (red lines). From top to 
bottom we have the profiles of gas temperature Tgas, density Pgas, 
entropy Sgaa = Tg^/nli^, (5 = (SksT) and surface 

brightness PgasTg^ . The two columns in each row are halos within 
different mass bins. The shaded regions represent the Icr scatter. 


ulations, we find that effective f{R) halos and ACDM halos 
have very close density and entropy profiles in this core region 
(see the second and third rows in Fig.|^. 

The shape of the gas density profile outside the core region 
(r > i?vir/2) can be described by the /3 model ll20l . where /3 
is the ratio between the specific kinetic energy /Ccdm (kinetic 
energy per unit mass) of cold dark matter and the specific in¬ 
ternal energy Sg^s (internal energy per unit gas mass) of gas 

. ^ Aicdm ^ ^ .2) 

~ ■^^gas -Pgas//0gas/(7 - 1) SfegT 

If we consider the 3-D velocity dispersion as a measure 
of the temperature of cold dark matter ED, then P is simply 
the ratio of temperatures between cold dark matter and gas. 

As shown in Fig. the temperature profile does not vary 
significantly outside the core region and so gas and cold dark 


matter can be roughly treated as isothermal. Furthermore, 
if we assume a hydrostatic equilibrium between them, the 
shapes of the gas and dark matter density profile are related 
byHH 


_ dlnpgas/dlnr 
dlnpcdm/dlnr 

In the ACDM model, this feature has been verified by a num¬ 
ber of numerical simulations 1231. As shown in the fourth 
row of Fig. 1^ it is interesting to find that /3 is almost identi¬ 
cal in f{R) and ACDM models. According to Eq. (|^, out¬ 
side the core region, the gas density profile in f{R) gravity 
should trace the dark matter density profile in the same way 
as that in ACDM. As indicated in the second row of Fig. 
at r > i?vir/2, in the standard catalog, the gas density pro¬ 
file in f{R) gravity matches the ACDM result (blue and black 
lines); in the effective catalog, the gas density profile in the 
f{R) model matches the ACDM one with a rescaled baryon 
fraction (red and green lines). 

According to the gas density and temperature profiles, it can 
be expected that the profiles of other physical quantities that 
depend on gas density and temperature should have the fol¬ 
lowing features: inside the core region, effective halos should 
behave similarly to ACDM halos; outside the core region, 
effective halos should resemble ACDM halos with rescaled 
gas fractions by /M^^\ For illustrative purposes, we 

show the profile of the surface brightness I{r). Since Hr) is 
related to Pgas(?’)^T'gas(r)i/2, in the last row of Fig.|^we only 
show the profile of Pgas(?’)^Tgas(?')^^^- It is evident that the 
Pgas(?')^Tgas(r')^/^ profile has the features as just described 
the above. 


V. GAS FIfACTIONS 


After studying the gas density profiles, we turn to the cu¬ 
mulative halo gas mass fraction. Since outside the core re¬ 
gions, the gas density profile in f{R) gravity traces the dark 
matter density profile, the cluster gas fraction is a tracer of 
the true, rather than the effective, halo mass. In objects as 
large as galaxy clusters, it is well known that the gas fraction 
is more sensitive to rather than the theory of gravity. 

This is illustrated in Fig.|^(also see Ref. Il24ll ). where we can 
see that the standard halos in f{R) gravity and ACDM ha¬ 
los have essentially identical gas fraction profiles (blue points 
versus black triangles). On the other hand, the gas fraction in 
the /(i?) effective halo agrees with that of ACDM halos with 
rescaled Of,. 


VI. SCALING RELATIONS 


We next turn to the scaling relations of cumulative physical 
quantities. The cumulative quantities, such as the luminosity 
Lx, describe the properties of clusters as a whole. In gen¬ 
eral, they are less sensitive to the profiles in the core region 
because the latter has a small volume and contributes little to 
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Standard Catalog Effective Catalog 



FIG. 3. Gas temperature as a function of halo mass in the ACDM model (black symbols) and f{R) gravity (red symbols). Left: the standard 
halo catalogue of f{R) gravity is used. Right: the effective halo catalogue of f{R) gravity is used. The shaded regions represent Icr scatter. 



logio (r/ [kpc/ h ]) logio (r/ [kpc/ h ]) 


FIG. 4. The cumulative halo gas fraction as a function of the distance from halo centre. In the standard f{R) halo catalog (blue points), the 
gas fraction agrees very well with that of the ACDM halos (black triangles). In contrast, the gas fraction of effective f{R) halos (red points) 
agrees with that of ACDM halos with rescaled Qb (green triangles). The left panel is for a mass bin in which most halos are unscreened and 
the right panel is for a mass bin in which halos are partially screened. 


the total cumulative value, but they are more sensitive to the 
profiles at a relatively larger radius. As noted earlier, outside 
the core region, the profiles of effective halos in f{R) grav¬ 
ity closely resemble those in the ACDM model with rescaled 
gas fractions. This feature, as we shall show later, can help us 
build a connection between the scaling relations of physical 
quantities of f{R) effective and ACDM halos. 


We start with the temperature-mass relation. As shown in 
the left panel of Fig. in the standard catalog, the gas tem¬ 
perature of unscreened f{R) halos are higher than those of the 
ACDM halos with the same masses, which is consistent with 
what was found in ll^ . However, in the effective catalog, as 
shown in the right panel of Fig. the Tgas-M relation in the 
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two models agree very well with each other 

. (4) 

This is a natural result from the gas temperature profile as 
discussed above. 

For gas densities, we have seen above that outside the core 
region we have 




^ _ ACDM 

»#/(«) 






where r^ore is the core radius. For an f{R) effective halo with 
the same mass as a ACDM halo, from the 

above equation and the fact that f{R) effective and ACDM 
halos have very similar temperature profiles, it follows that 



where a and b are indices of power. 

The above relation is one of the key results in this paper and 
it indicates the important relation of cumulative gas quantities 
between f{R) effective halos and ACDM halos. In order to 
test the validity and illustrate the use of Eq. Q, we investigate 
three most important and frequently used quantities in cluster 
surveys as examples. The first one is the x-ray luminosity Lx 
which, for a cluster, can be written as 


dr47rr^(p^ 


1 ^ a ^ b 




Ax(< r)= [ drATTr"^. (6) 

Jo 

Lx is sensitive to the details of the gas distribution in the cen¬ 
tral region and depends on the dynamical state of the cluster. 
As a cluster mass proxy, Ly^-M has large scatters. 

The second one is the integrated SZ Compton y-parameter 

Ysz{< r) = [ drAirr'^Pe , (7) 

meC‘^ Jo 

where ar is Stefan-Boltzmann constant, me the electron mass 
and c the speed of light. Pe is the electron pressure, which is 
given by Pe = ^ngasA:BTgas • 

The third one is the x-ray equivalent of the integrated SZ 
flux, the Yx parameter IfTTlI . 


^x(< r) = Tgas [ dr47rr^/^gas , (8) 

where Tgas is the average mass-weighted temperature. Ix- 
M relation is relatively insensitive to the dynamical state of 
clusters and to the detailed modeling of gas physics and 
the Yx-M scaling relation is often practically used in x-ray 
surveys (e.g. ESI). 


TABLE I. The average relative differences of scaling relations be¬ 
tween f{R) and ACDM models. 



Before Rescaling 

After Rescaling 

< L«"VlACdm > 

27.2% 

13.1% 

< > -1 
< y^(fi)/yACDM ^ 

19.7% 

19.6% 

3.1% 

3.2% 


From Eq. Q, it follows that 

« yACDM(^ACDM ^ ^ 

(9) 

The numerical results of the scaling relations are shown in 
Fig. [5] from which we can see that without the rescaling of 
Eqs. (j^, the f{R) and ACDM models have very different 
scaling relations (left panels), but after the rescaling not only 
the mean value but also the scatters of the scaling relations 
in the two models highly resemble each other (right panels). 
In particular, for Ysz and lx after the rescaling the average 
relative difference between the models are only at a level of 
3% (see Table|I|. Therefore, using this rescaling method, the 
complicated effect of modified gravity on the cluster scaling 
relations can be accurately modelled and the error is much 
smaller than the typical uncertainty caused by the modelling 
of baryonic physics in galaxy formation (e.g. ||6l). 


Mf(R) ^sz ) 


M, 


/(fl) 


V 

M/(^) J 

rf{R) 




M] 


Eff tz/(-R) 




VII. DISCUSSION 

In this paper, we investigated the gas properties in a rep¬ 
resentative non-standard gravity model, f{R) gravity, based 
on a suite of non-radiative hydrodynamical simulations. We 
studied both the profiles and the cluster scaling relations of 
gas properties. We found that the effective halos, which are 
proposed in ifTSl . in f{R) gravity have similar temperature 
profiles as ACDM halos with the same masses. For the gas 
density profile, effective halos closely resemble ACDM halos 
in the core region, while outside the core region they behave 
like ACDM halos with rescaled gas fractions. 

Based on those observations, we have demonstrated that 
not only the mean value but also the scatters of the scaling 
relations of cumulative gas quantities in f{R) effective halos 
with rescaled gas mass fractions are very similar to those in 
the ACDM model, cf Eq. 0. In particular, for Fsz and Yx, 
which are the two most frequently used mass proxies in cluster 
surveys, our rescaling method enables us to model the effect 
of modified gravity on the scaling relations with an accuracy 
of ^ 3%, which is much smaller than the typical uncertainty 
caused by the modelling of baryonic physics in galaxy forma¬ 
tion (e.g. 161). The error in modelling the impact of modihed 
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FIG. 5. Cluster x-ray luminosity (upper panels), Ysz parameter (middle panels) and lx parameter (lower panels) as a function of halo mass 
in the ACDM model (black symbols) and f(R) gravity (red symbols). Left: results without the rescaling of Eqs. ij^. Right: results with the 
rescaling of Eqs. 12}. The results are almost identical in the two models in the latter case while differences are apparent in the former case. 
The shaded regions represent Icr scatter. 


gravity on the scaling relations is a minor source of the overall 
uncertainty in cluster cosmology. This, therefore, provides an 
accurate way to calibrate the scaling relations in f{R) gravity 
and a way to correctly interpret them in the context of non¬ 
standard gravity theories. It will enable us to test gravity using 
cluster observations in a fully self-consistent way, and avoid 
substantial systematic biases. Moreover, this provides a reli¬ 
able way to analyze the gas physics in modified gravity using 
much more efficient pure cold dark matter simulations with 
the known knowledge in the ACDM model. 

Although we illustrated our idea using a specific f{R) 
model with a fixed value of fuQ, our main results are ex¬ 
pected to apply to other values of fuQ and other f{R) models 
as well. This is because the f{R) model with the value of fno 
chose here contains the totally screened, partially screened, 
and totally unscreened halos. As shown in Fig.|^ our rescal¬ 
ing method works very well for all these halos that have dif¬ 
ferent levels of screening. The difference in different models 
of /(i?) gravity and different values of fj^Q only lies in the 
exact mass ranges within which halos are totally unscreened. 


totally screened or partially screened. Therefore, varying fuQ 
does not affect the overall workability of our method and our 
method works for other f{R) models as well. 

From the above argument, our main results are also ex¬ 
pected to apply to other models that employ the chameleon 
screening mechanism EH as well as dilaton Il28l . sym- 
metron ll29l models and their generalisations OOl . Therefore, 
our method can have much wider applications and provide a 
useful way to constrain modified gravity theories using obser¬ 
vations from upcoming cluster surveys such as eROSITA 0. 

Acknowledgments J.H.H. acknowledges support of 
the Italian Space Agency (ASI), via contract agreement 
1/023/12/0. B.L. acknowledges support by the UK STFC 
Consolidated Grant ST/L00075X/1 and RF040335. This 
work has used the DiRAC Data Centric system at Durham 
University, operated by the Institute for Computational 
Cosmology on behalf of the STFC DiRAC HPC Facility 
(www.dirac.ac.uk). This equipment was funded by BIS 
National E-infrastructure capital grant ST/K00042X/1, STFC 
capital grant ST/H008519/1, STFC DiRAC Operations grant 









7 


ST/K003267/1 and Durham University. DiRAC is part of the National E-Infrastructure. For access to the simulation data 

and halo catalogs, please contact J.H.H. 


[1] A. Albrecht, et. al., arXiv:astro-ph/0609591 

[2] R. A. Sunyaev, and I. B. Zel’dovich, Annu. Rev. Astron. Astro- 
phys., 18, 537 (1980). 

[3] A. Merloni, eROSITA Science Book, arXiv:1209.3114 

[4] R. Laureijs, Euclid Definition Study Report, arXiv:1110.3193 

[5] S. Ameglio, S. Borgani, E. Pierpaoli, K. Dolag, S. Ettori, and 

A. Morandi, Mon. Not. R. Astron. Soc., 394, 479 (2009). 

[6] D. Eabjan, S. Borgani, E. Rasia, A. Bonafede, K. Dolag, G. Mu- 
rante, L. Tornatore Mon. Not. R. Astron. Soc., 416, 801 (2011). 

[7] C. Loken, M. L. Norman, E. Nelson, J. Burns, G. L. Bryan, 
P. Motl, Astrophys. J., 579, 571 (2002); S. Borgani et. al, 
Mon. Not. R. Astron. Soc., 348, 1078 (2004); M. Roncarelli, 
S. Ettori, K. Dolag, L. Moscardini, S. Borgani, G. Murante, 
Mon. Not. R. Astron. Soc., 373, 1339 (2006); J. H. Croston, 
et. al., Astron. Astrophys., 487, 431 (2008). 

[8] E. Schmidt, A. Vikhlinin, W. Hu, Phys. Rev. D80, 083505 
(2009); M. Cataneo, D. Rapetti, F. Schmidt, A. B. Mantz, 
S. W. Allen, D. E. Applegate, P. L. Kelly, A. vonderLinden, 
and R. G. Morris, Phys. Rev. D92, 044009 (2015); A. Barreira, 

B. Li, W. A. Hellwing, L. Lombriser, C. M. Baugh, and S. Pas- 
coli, J. Cosmo. Astropart. Phys., 04, 029 (2014). 

[9] W. Hu, and I. Sawicki, Phys. Rev. D76, 064004 (2007). 

[10] J. -h. He, B. Li, A. J. Hawken, B. R. Granett, Phys. Rev. D90, 
103505 (2014). 

[11] A. V. Kravtsov, A. Vikhlinin, D. Nagai, Astrophys. J., 650, 128 
(2006). 

[12] B. Li, G. -B. Zhao, R. Teyssier, and K. Koyama, J. Cosmo. As¬ 
tropart. Phys., 1, 51 (2012). 

[13] R. Teyssier, Astron. & Astrophys., 385, 337 (2002). 

[14] J. -h. He, and B. Wang, Phys. Rev. D87, 023508 (2013). 

[15] P. Ade et al., Astron. Astrophys., 571, A16 (2014). 

[16] E. Bertschinger, arXiv:astro-ph/9506070, S. Prunet, C. Pichon, 
D. Aubert, D. Pogosyan, R. Teyssier, and S. Gottloeber, Astro¬ 
phys. J. Suppl., 178, 179 (2008). 


[17] S. R. Knollmann, and A. Knebe, Astrophys. J. Suppl., 182, 
608 (2009); S. P. D. Gill, A. Knebe, and B. K. Gibson, 
Mon. Not. R. Astron. Soc., 351, 399 (2004). 

[18] J. -h. He, A. J. Hawken, B. Li, and L. Guzzo, Phys. Rev. Lett., 
115, 071306 (2015). 

[19] J. -h. He, B. Li, and A. J. Hawken, Phys Rev. D92, 103508 
(2015). 

[20] E. Komatsu, and U. Seljak, Mon. Not. R. Astron. Soc., 327, 
1353 (2001); Y. Suto, S. Sasaki, and N. Makino, Astron. J., 
509, 544(1998). 

[21] E. Komatsu, and U. Seljak, Mon. Not. R. Astron. Soc., 336, 
1256 (2002). 

[22] H. Mo, F. Van den Bosch, and S. White, Chapter 8 in "galaxy 
formation and evolution", Cambridge University Press, New 
York, (2010). 

[23] J. F. Navarro, C. S. Frenk, and S. D. M. White, 
Mon. Not. R. Astron. Soc., 275, 720 (1995); G. F. Lewis, 
A. Babul, N. Katz, T. Quinn, L. Hemquist, and D. H. Wein¬ 
berg, Astrophys. J., 536, 623 (2000); K. Yoshikawa, Y. P. Jing, 
and Y. Suto, Astrophys. J., 535, 593 (2000). 

[24] B. Li, J. -h. He, and L. Gao, Mon. Not. R. Astron. Soc., 456, 
146 (2016). 

[25] C. Arnold, E. Puchwein, and V. Springel, Mon. Not. R. As¬ 
tron. Soc., 440, 833 (2014). 

[26] A. Vikhlinin, R. A. Burenin, H. Ebeling, and W. R. Forman 
et. al, Astrophys. J., 692, 1033 (2009). 

[27] J. Khoury, and A. Weltman, Phys. Rev. D69, 044026 (2004); 
J. Khoury and A. Weltman, Phys. Rev. Lett., 93, 171104 (2004). 

[28] P. Brax, C. van de Bruck, A. -C. Davis, and D. Shaw, 
Phys Rev. D, 82, 063519 (2010). 

[29] K. Hinterbichler, J. Khoury, Phys. Rev. Lett., 104, 231301 

( 2010 ). 

[30] P. Brax, A. -C. Davis, and B. Li, Physics Letters B, 715, 
38 (2012); P. Brax, A. -C. Davis, B. Li, and H. A. Winther, 
Phys. Rev. D86, 044015 (2012). 



